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A numerical method for the calculation of electronic structure of a nanosystem composed of a 
pseudoisocyanine (PIC) molecule assembled on a silver nanoparticle is developed. The electronic 
structure of the silver nanoparticle containing 125 atoms is calculated within the local density 
version of the density functional method. A model of an Ag atom embedded in the center of a 
spherical jellium cluster is used. The host electron Green's function is calculated by means of 
the spherically symmetric expansion. The principal theoretical tool is the scattering theory using 
the Green's-function method. The molecule - silver nanosystem interaction is studied using the 
approach similar to that of the Anderson model for transition metal impurities in solids. Localized 
levels are shown to split off from the top of the band of the Ag nanoparticle. The electronic structure 
calculations yield information on the character of chemical bonding in the PIC molecule — silver 
particle nanosystem. 



PACS numbers: 

I. INTRODUCTION 

Organic-inorganic nanohybrid materials that utilize 
noble metals (silver or gold) , functionalized with organic 
or biological constituents can produce unique physical 
and chemical properties that otherwise are not possi- 
ble in single component materials** These materials have 
given rise to growing interest in theoretical methods that 
can calculate the electronic structure and transport prop- 
erties of nanoscale devices A major problem with 
these nanostructures is to optically control molecular self- 
organization on metallic surface. J-aggregates of cyanine 
dyes make a fascinating topic of research due to their 
outstanding optical properties^ 5 Spectroscopic peculiari- 
ties are the result of exceptionally strong electronic in- 
teractions between the transition dipole moments of the 
dyes that give rise to extended exciton states after photo 
excitation^ The excitonic optical spectrum depends on 
the details of the structural and electronic arrangement 
of the dye molecule^ 

Adsorption of dyes to nanoparticles of noble metals 
presents a special interest^ Study of the interaction of 
an adsorbed dye, in both its ground and excited states, 
with the energy states of the conduction band of metal 
is also of interest and may rely on theoretical calcula- 
tions of the energy electronic structure of the nanosystem 
made of pseudoisocyanine (PIC) assembled on a silver 
nanoparticle. Moreover, the nature of coupling between 
excitonic molecular J-aggregates and a metallic nanopar- 
ticle is still not completely understood, since to the best 



of our knowledge no attempts had been made to calcu- 
late the electronic structure of Ag + PIC nanosystems 
from the "first principles" In addition, a major prob- 
lem of molecular engineering is to control molecular self- 
organization. 8 This process is governed by the formation 
of chemical bonds, e.g. covalent or hydrogen bonds. The 
assembly process on nanoparticles is also affected by the 
molecule — substrate interactions. Therefore, under- 
standing bonding between the molecules and nanopar- 
ticle is crucial in order to be able to choose appropriately 
the molecular and substrate materials for a nanosystem 
design. On the theoretical side, the new computational 
techniques allow one to predict the preferential geometry 
of the nanostructure arrangements as well as the strength 
and chemical bonds involved directly from the fundamen- 
tal quantum mechanical laws. 

We demonstrate how ah initio calculations allow one 
to explain the shape of the observed superstructures, to 
elucidate the role of electronic structure and the molecule 
- silver particle bonding and to reveal details of the 
nanosystems not yet experimentally accessible. 

Quantum electronic-structure calculations allow us to 
understand the macroscopic properties of complex poly- 
atomic systems (specifically organic PIC molecule assem- 
bled on a silver nanoparticle) in terms of the microscopic 
states available to the electrons described by their wave 
functions and the nanosystem electron density p nano (r). 
For these larger systems it is frequently of a greater inter- 
est to know the change in the electronic structure associ- 
ated with the changes in the Kohn-Sham 9 effective poten- 
tial V e f f (r) and the electron density. Our work develops 
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a theoretical technique that makes the direct calculation 
of such changes possible. This technique is based on the 
self-consistent Green's function method^ and standard 
density-functional theory (DFT) in combination with the 
local- density approximation (LDA)^ which maps the 
many-electron interaction problem onto a self-consistent 
treatment of non-interacting quasi-particles moving in an 
effective potential 

/ nano /_/\ 
| r _r' dr/ + ^ c(r) ' (1) 

where V ex t(r) is the external potential and V xc (r) is the 
local exchange-correlation potential. 

II. FORMULATION OF THE PROBLEM 

We present in this paper results of numerical calcu- 
lations of the electronic structure of nanosystem: PIC 
molecule + small Ag particle, whose electronic structure 
is described by the Hamiltonian 

H = HjeiUum + Hid + H PIG + HpIC-Ag + ^PIC-Ag (2) 

The first term stands for the spherical jellium cluster 
Hamiltonian (see below) made of 5s-valence electrons of 
Ag-atom. The next term in the Hamiltonian describes 
the subsystem of Ag 4e?-electrons embedded in the spher- 
ical jellium cluster. Hpjc stands for the PIC molecule 
Hamiltonian. The last two terms of the Hamiltonian @ 
describe the resonant(r) (due to the presence of the Ad 
state of Ag atoms) and potential (p) (due to the Ag 5s 
states) scattering. This distinction between the two types 
of scatterring by atomic potential is discussed in detail 
in Ref. [H (see also[l(J. 

The method is based on the separation of the Hamil- 
tonian ((2|) into the two parts 

H = H° LDA + AH (3) 

Here Hl DA = Hj e m um + H^d is the Hamiltonian of 
the reference system, which includes an effective single- 
particle LDA potential V° DA (r). AH = H PIC + 
Hpic-A g + Hpic~Ag * s the Hamiltonian of the "differ- 
ence" system described by the potential AV(r), which is 
the difference between the self-consistent Kohn-Sham (fT]) 
effective potential V e //(r) and V^ DA (r). 

Four different computation schemes are used, each 
of which providing mutually complementary informa- 
tion about the electronic spectra of complex nanosys- 
tem. The first one is the spherical jellium cluster model 
with the central Ad Ag atom embedded in the center of 
a silver nanoparticle-^ The second method is based on 
the Green's function approach. Then the physical ef- 
fect shows up exclusively in variations of the Kohn-Sham 
effective potential V e //(r)2 and of the nanosystem elec- 
tron density p nano (r). In this approach, we follow a two- 
step concept and split the entire problem into two parts, 



where each of them is far less complicated than the orig- 
inal onei 14 ' 15 The whole nanosystem is decomposed into 
the reference system characterized by the ground-state 
electron density P°ld A {y) and the difference system with 
the electron density 

Ap(v) = p— (r)-p° LDA (r), (4) 

which is the difference between the self-consistent elec- 
tron density p nano (r) of the nanosystem and the ground- 
state electron density P lda( t ) corresponding to the 
Hamiltonian 

Hl DA = -V 2 + Vi DA (v) (5) 

(atomic Rydberg units are used). Here H^ DA = 
H jellium + Hid- The complete set of electronic states 
of the reference system is represented by the reference 
Green's operator 

G°(z) = (z-H° LDA )- 1 (6) 

depending on the complex variable z = e + in (r/ > 0). 
The calculation of G°(z) along a properly chosen contour 
in the complex energy plane defines step one of our ap- 
proach. As the second step, the Green's operator G(z) of 
the nanocluster-assembled system (and the charge den- 
sity related to it) is determined by solving the Dyson 
equation 

G(z) = G°{z) + G°(z)AVG{z) (7) 

in a self-consistent way. A detailed description of our 
approach to solving Eq. (J7J is presented in Appendix. 

The third semiempirical method ZINDO/Si^ was used 
for the calculation of electronic structure of a neutral PIC 
molecule. The fourth DFT/B3LYP/6-31G^ method was 
applied to the equilibrium geometry calculation of PIC 
molecule by the standard DFT-method. 

We would like also to indicate here the limitations of 
our approach to the treatment of Ag + PIC nanosystem 
based on the spherical jellium model for the Ag nanopar- 
ticle. The previous structural researched using method 
of the molecular dynamics shows that very diverse non- 
spherical geometries of the nanoparticle are possible. Re- 
laxation of the PIC molecule and nanoparticle recon- 
struction in response to the PIC molecule are also not 
taken into account. Including all the processes into the 
calculation scheme can strongly complicate the procedure 
and we leave it for our future research. 



III. CALCULATION OF THE ELECTRONIC 
STRUCTURE OF AG - NANOPARTICLE 

First, we consider the electronic structure of silver 
nanoparticle. Cluster calculations are traditionally em- 
ployed in studies of surface and bulk materials. They 
help us to understand how the physical properties evolve 
from a free atom to a finite-size system. In recent years, 
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microclusters on the basis of Ag atoms attracted a lot 
of interest, due to the growing technological significance 
of nanosystems including Ag and Au nanoparticles . 19 ' 20 
A variety of theoretical models have been proposed for 
the calculations of cluster electronic structure. However, 
the most precise models allow one to obtain their elec- 
tronic structure with a small number of atoms only. For 
large atomic aggregations, these methods cannot be ap- 
plied successfully, and simpler models are employed in 
this case. Thus, for the study of electronic properties of 
sp-bonded metal clusters, a jellium model is used. 22 But 
the models based on the jellium approximation are not 
directly suitable for the investigations of metal clusters 
containing atoms with localized d shells. As a conse- 
quence a model of an Ag atom embedded in the center 
of a spherical jellium cluster is applied for a description 
of a small Ag metallic nanoparticle containing localized 
Ad electrons. The DFT approach in the local density ap- 
proximation was used in computation! 13 ' 21 Each atom of 
the nanoparticle is mimicked by a single atom embedded 
in the center of a jellium sphere with the ru radius, de- 
termined according to the position of the real atom with 
respect to the cluster surface, is then the shortcut 
between the atom and the cluster surface. 

The electronic structure of an Ag atom embedded in 
the jellium sphere is obtained within the DFT framework 
from a self-consistent solution of the Kohn-Sham equa- 
tions (in atomic Rydberg units)£i2i 



VE D A(r)]ipm(r) 



(8) 



where 
V£ DA (r) 



with the Vosko et al forrn^ 3 - of the local exchange- 
correlation potential V xc (r). Here e n i and ip n i are single 
electron energies and wavefunctions, respectively; Z° is 
the nuclear charge of the Ag atom. The electron density 
of the jellium cluster with an Ag atom is 
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"( r ) =^2fnl I -0ni(r) 



(10) 



where the coefficients /„; are the occupation numbers of 
the states with quantum numbers n, I, and the summa- 
tion is over all states of the atom-in- jellium nanoparticle. 
The radial distribution of the positive jellium background 
is given by 

p+(r) = [3N val (N at - l)/47rr| f ]6(r M - r), (11) 

where 0(x) is the unit step function, N at is the number 
of atoms in the cluster (including the specific Ag atom) , 
N va i is the number of valence electrons in the Ag atom. 
The cluster radius rjvf is found from the expression 



where r c is the Wigner-Seitz radius. We have used here 
r c = 3.02356 a.u., N at = 125, and N val = 1. The nu- 
merical integration of the Kohn-Sham equation for the 
Ag atom in a jellium sphere is carried out by means of 
the Milne method. 2 - The free Ag atom eigen- energies e n i 
and wave-functions ip„i are calculated by means of the 
semirelativistic RATOM program^ The self-consistency 
procedure for the potential V® DA (r) is carried out in a 
mixed fashion. The first two iterations use the arithmetic 
average scheme, which later on is effectively substituted 
by the Aitken scheme. 26 

The electronic ground-state configurations of the jel- 
lium sphere containing the central Ag atom were deter- 
mined in the following way. The rules for the energy 
level occupation separately in jellium and in a free atom 
are well known. Clearly, insertion of an atom into the 
jellium sphere center does not change the atom and jel- 
lium field symmetry, and, consequently, the symmetry 
of their electronic states is not changed either. One can 
suppose therefore that the number of electron states with 
the same symmetry in the jellium sphere with the central 
atom is equal to the sum of such symmetry states of the 
jellium and the atom. The sequence of energy levels of 
atom-in-jellium is obtained by solving the self-consistent 
equations (J5J - (|10p for the different angular I and princi- 
pal n quantum numbers. The energy levels are occupied 
in accordance with the Pauli principle. The highest levels 
can be partially occupied. 

We have computed the total energy of the jellium 
spheres with the Ag atom in the center for various occu- 
pation numbers of the upper levels. For the calculation 
of the total energy we employed the equation 
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Pout( r ) X 



p~ u t (r') - 2pr (r') 



dr' 



dr. (13) 



r M = Nl( 3 r, 



(12) 



Here the indices in and out indicate the input and output 
data of the latest self-consistency iteration, respectively; 
e°"*(r) is the exchange correlation energy density of a 
homogeneous electron gas with the density Po Ut (r) pa- 
rameterized according to Vosko et ali& 

The calculations have shown that the spherical jellium 
nanoparticle with the central Ad Ag atom has the follow- 
ing energy spectrum. The eight lowest energy levels are 
identical to those of the Ad atom core; the rest of them 
are similar to the states in a spherical potential well. In 
this paper the quantity n = n n + 1 has been considered 
as the principal quantum number; n n is the number of 
nodes of the wave-function of the corresponding energy 
level. Thus, the electronic configuration of the central Ag 
atom in the jellium sphere is ls 2 2s 2 lp 6 3s 2 2p 6 ld 10 4s 2 3p 6 . 
The electronic states of the 125 atoms are arranged in 
the following order: 2d 10 l/ 14 l 5 18 l/i 22 2/ 14 li 26 lj 30 . The 
self-consistent potential V® DA (r) (Eq.@) with the elec- 
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FIG. 1: (Color online) A self-consistent Ag-nanoparticle po- 
tential of the jellium sphere containing the central Ag atom 
(the radial distribution, in atomic units). The states of sub- 
shell 2d 10 are split into 4d 5/2 and 4d 3/2 levels due to the spin- 
orbit interaction. 



tronic levels for an Ag atom embedded in a jellium sphere 
is presented in Fig. Q] 

An important characteristic of the Ag nanoparticles is 
their ionization potential (IP). It has been shown in Ref. 
l27l that a strong correlation exists between the chemisorp- 
tion reactivity of a small transition metal cluster and its 
ionization threshold. Here we report the IP's for the Ag 
nanoparticles containing up to 160 atoms. In these cal- 
culations we simulate atoms of an Ag nanoparticle by Ag 
atoms embedded in the center of jellium spheres of vari- 
ous sizes. The jellium sphere radius tm is denned by the 
short cut between the atom and the cluster surface. 

The IP of an atom- in- jellium was obtained using the 
ground state theory (LDA-method) by self-consistent cal- 
culation according to Eq. (fTBf and subtraction of the 
total energies of neutral and ionic ground states. It was 
found that the ionization thresholds of atoms of the Ag 
nanoparticles differ and depend on the position of the 
atoms with respect to the cluster surface. The lowest 
atomic ionization threshold has been chosen as the IP 
of the Ag nanoparticle. For the neutral Ag atom we 
obtained IP^ = 7.96eV and IP^ m = 7.576eV. In- 
stead of calculating the IP's from the data as the total 
energy differences between the neutral cluster and ions, 
we use an alternative method, based on the Slater transi- 
tion state approach^ and Janak theorem^ for the den- 
sity functional theory, which allow us to calculate the 
excitation energy of adding (removing) an electron to 
(from) the system from (to) the infinity. This scheme 
can be derived in the following way. The total energy 
E to t difference between the final and initial states for the 
process of electron addition to the one-electron state j 
can be calculated as an integral of total energy deriva- 
tive with respect to the occupancy qj. This derivative, 



Ej = dE/dqj, corresponds to the Kohn-Sham eigenvalue 

1 

E( Qj = 1) - E( qj = 0) - J dq^qj) ~ £, (0.5). (14) 



Eq. (Thij) becomes exact if the LDA eigenvalue £j(qj) is 
a linear function of the occupancy, which is usually true 
to within a good accuracy. 

The vertical (or adiabatic) ionization potentials^ of 
the clusters are evaluated from the highest occupied 
molecular orbital (HOMO) energy of the neutral clus- 
ters. The calculated IP of an Agi25 nanoparticle is 
Ip theor = 4 73e v. When the number of atoms in the 
nanoparticle increases the IP converges to the work func- 
tion $ of the corresponding metallic half-space. The 
work function of 4.46 eV was determined experimentally 
for Ag(lll)4i Our calculation yields the value of 4.73 
eV for IP t 1 h 2 e 5 or . The LDA HOMO-LUMO (lowest un- 
occupied molecular orbital) gap between occupied (e^o) 
and unoccupied (s g is) orbitals for the Ag nanoparticle is 

A HOMO-LUMO = £u2(> _ ^ = X 58eV 

It is well-known that Ag aggregates are formed in a 
colloidal aqueous solution of NaCljl Therefore studying 
electronic properties of small metallic particles in differ- 
ent dielectric matrice s 32 ' 33 ' 44 presents a great scientific 
and practical interest as well as those of free Ag nanopar- 
ticles. When a metallic particle is placed in a dielectric 
medium, polarization changes p po i (r) are induced on the 
particle surface, which produce the potential 

Vpoi{r)=2jf^±dv> (15) 

In Ref. HU, we have shown that the ground state of a 
metal nanoparticle embedded in a dielectric matrix can 
be described by the self-consistent Kohn-Sham equations 
with the effective potential 

V s (v) = V(r) + V pol (r) (16) 

where V(r) has the same form as the effective potential 
for the nanoparticle in vacuum. Averaging V po i(r) and 
substituting the result into (fTB"]) yields 

f V(r) + ^V(r M ) r < r M , 
V e {r)={ 

{ V(r)/e r > r M . 

The calculations for the Ag\2^ nanoparticle embedded 
in a dielectric medium with the relative permeability 
e = 61.1 (50% water solution of NaCl) have shown that 
with the increasing dielectric permeability e the poten- 
tial profile near the jellium edge becomes steeper. The 
bottom of the potential well Vbottom and the single elec- 
tron energy levels Ei rise when the cluster is embedded 
in a dielectric matrix. It should be emphasized that the 
oscillations^ 3 - in the electronic density of jellium cluster 
in vacuum are suppressed in the dielectric media. The 



5 



amount of the electronic charge beyond the jellium edge 
(electronic "spill out") increases with the increasing e. 
This is caused by a positive shift of the electronic eigenen- 
ergies of Ag nanoparticles in a dielectric medium and an 
extension of the corresponding wave functions. 

The static dipole polarizability 3 -^ is also well-known 
to be intimately related to the electronic structure of a 
nanoclusteri^ We discuss here the photo-response of an 
isolated Ag nanoparticle to an external electromagnetic 
field in terms of the frequency-dependent polarizability 35 

oo 

a(u) = -y J dr'r' 3 6p(r',Lu) (17) 
o 

where Sp(r' ,u>) is the change in the charge density. The 
calculation of the static dipole polarizability a(uj = 0) 
for an Ag nanoparticle is carried out within the time- 
dependent local-density approximation (TDLDA) 3 ^ us- 
ing the self-consistent solution of the set of equations 

5p(r,u) = f X o(r,r';u)SV(r>)dr', (18) 



SV(r,u>) = V ext (r,uj) + V ind (r,u), (19) 

and 

^ (r ,„) = 2 /«^|dr' + ^£)«r,„). ,20) 

Here V ext (r,uj) and V in d{r,uj) are the external field 
with the frequency uj and the induced field, respec- 
tively; Xo( r , r '; w ) is the susceptibility function in the 
independent-particle approximation. For tm — 15.1178 
a.u. (Agi2 5 cluster) a(uj = 0) = 4.82A 3 / atom. The bulk 
atomic polarizability is 4.33 A 3 /atomj2£ 

IV. RESULTS AND DISCUSSION 

We present here a theoretical approach to the organic 
molecules interacting with silver in terms of numerically 
solvable DFT-models. In practice, these models apply 
to nanosystems (PIC-molecule on Ag-particle) and pro- 
vide an understanding of doping silver nanoparticles by 
PIC adsorbates. The chemical structure of PIC molecule 
is exhibited in Fig. [3] The density-functional calcula- 
tions are performed in order to determine the equilibrium 
structure of these PIC molecules by the DFT/B3LYP/6- 
31G-methodii£ The molecule geometries are optimized 
using the B3LYP-exchange-correlation energy functional 
and potential. The minimum basis set is capable of pro- 
ducing reasonable results is 6-31G. Since the geometry 
of PIC within the Ag+PIC-nanosystem is undoubtedly 
disturbed as compared to the PIC monomer and possibly 
resembles a PIC aggregate, these calculations were per- 
formed as follows. First, the equilibrium geometry of the 




FIG. 2: (Color online) Configuration of the PIC molecule 
calculated by DFT/B3LYP/6-31G-method: H(black) stand 
for the hydrogen atoms; C(red) stands for the carbon atoms; 
N(yellow) stands for the nitrogen atoms; F(blue) stands for 
the fluorine atom. 



tetramer (PIC:C1)4 was found. After that, the averaged 
geometry of the monomer (PIC:C1) was obtained by av- 
eraging over four molecules of the tetramer. The search 
of the equilibrium geometry of tetramer was performed 
by standard density functional method DFT/B3LYP/6- 
31G using GAMESS program packaged The electronic 
structure of a neutral PIC molecule was calculated by 
ZINDO/S-method^ 

It is assumed that the PIC molecule is adsorbed on 
the spherical surface of an Ag nanoparticle, setting down 
over the sphere so that its central carbon atom (the 
molecule center of gravity) is positioned in the point 
(2.223,0.680,0.705)Awith respect to the center of Ag 
nanoparticle. In the beginning we present the results 
of calculation of the electronic structure of an Agi 2 5 
nanoparticle obtained by the model of an atom in the 
center of a spherical jellium sphere and DFT in the local 
density approximation. Fig. [T] presents the energy level 
structure of the Ag jellium nanoparticle with 134 valence 
electrons (the Ag atom core electrons orbitals lie below 
e 3p 3/2 = - 4.27 Ry ). 

Fig[3] shows the Ag bound states by the red line. The 
peak heights are proportional to the degeneracy. The 
continuum density of states An(e) is given for the pos- 
itive energies corresponding to the delocalized (unoccu- 
pied) states by Eq. (|A18|) . Fig. [3] exhibits the contribu- 
tion of potential scattering [Eq. ([A17l) ] to the formation 
of Ag related bound states (e2g =-1.34 eV) and fall in the 
broad range A HOMO ~ LUMO . This level is obtained from 
Ea. (|AT7| with M substituted for Q- 1 [see Eq. (|ATT|) ]. 
The states 2d 10 (4d 3 / 2 + 4d 5 / 2 )l/ 14 l 5 18 l/i 22 2/ 14 H 26 lj 30 
in the occupied part of the spectrum are the Ag reso- 
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nances in the band (red line in Fig.Q and lie in the in- 
terval (-10,-5) eV. The resonant £29 level [see Figj3] arises 
at the energy 0.24 eV below the LUMO-level. Both the 
Ag resonances and E2 a states are found in the solution of 
Eq. (|A17l) with the full self-energy M(z). 

The calculation of Agi25 particle embedded in a dielec- 
tric medium with relative dielectric permeability e = 61.1 
(50% water solution of NaCl) have shown that the elec- 
tronic "spill out" increases in the dielectric medium as 
compared to vacuum and is equal for this case to 21.6 
electrons. 

The 26.14 eV width band of Ag+PIC nanostructure is 
formed by strongly hybridized C(s,p) (blue line), N(s,p) 
(yellow line), H(s) (black line) and Ag (red line) states. 
The unoccupied states lie above e LUMO = —1.58 eV. 
Furthermore, we find the acceptor-like £i g states almost 
entirely localized in the adsorbate PIC molecule and hy- 
bridized with the hydrogen Is levels. As discussed in 
Ref. Hfl the interaction between the adsorbate and the 
transition metal surface can be described as a two-state 
problem (adsorbate state and the d-band) leading to for- 
mation of bonding and anti-bonding states. Thus, an up- 
shift of the d-states should increase the adsorbate-metal 
interaction, since it would lead to the formation of an 
anti-bonding orbital closer to the Fermi level. Strong 
features appearing between -5 and -10 eV represent for- 
mation of a bonding orbital through the interaction of 
hydrogen Is state with the metal d-band, and this for- 
mation is typical for all transition metals. The bonding 
Ag 2d- Is states lie around -9 eV. The anti-bonding or- 
bital lies around e LUMO level and is hybridized with the 
resonance S2 g states. 

We will interpret the density of states (see figure [3]) 
of nanosystem Ag+PIC species basing on their bond- 
ing properties obtained from the Green's function cal- 
culations and simple molecular orbital considerations. 



For a more detailed study of the chemical bonding in 
these nanosystem we have calculated the electron-density 
change Ap(r) (see figure QJ. 

The electronic structure obtained for the Ag+PIC 
nanosystem has a terminal Ag — C bond. Similarly to all 
transition metals 4 ^ the C — C bond is known to be much 
stronger than the Ag — C bond. Thus, it is not energeti- 
cally favorable for Ag to enter the C chain, since this will 
break the stronger C — C bonds and form weaker Ag — C 
bonds. The Green's function calculations indicate that 
the excited Ag atoms in the nanosystem are involved in 
the chemical bonding. Simple molecular orbital consider- 
ations suggest that the C atom directly bonded to an Ag 
atom must be in an sp-hybridized state in order to more 
effectively interact with the Ag particle orbitals. And the 
Ag 2d orbitals can interact with the other p orbitals of 
the C atoms to form additional bonds. Thus, our cal- 
culations show that the Ag(d,s) — C(s,p) bond is of a 
considerable importance for understanding the nature of 
chemical bonding of Ag+PIC nanosystem. 

As shown in FigJ3J all Ag states are situated be- 
tween the C-states. From An(e) distribution as well 
as from its partial components An s (e), An p (e) and 
And(e), one can conclude, that d — s,p resonance is 
observed in (Ag+PIC)-nanosystem, was revealed in the 
compounds of metals, containing the filled d-shells, 
with non-metals.— In this case,the Ag 5s-type states 
(l/ 14 l 5 18 l/i 22 2/ 14 li 26 lj 30 ) form cr bonding and a* anti- 
bonding orbitals with the C sp- hybridized orbitals. At 
the same time, the 2d7r orbitals in Ag can interact with 
the C — C 7r orbitals. The resonance of the HOMO e^g 
with the hydrogen ls-level (cr*-symmetry) is the prin- 
ciple interaction in the Ag5s — His bond. This bond 
plays the crucial part in determining the chemisorption 
reactivity^ 

There are two factors that affect this orbital interac- 
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tion, first the HOMO energy, and second and more im- 
portantly, the symmetry of the HOMO. It follows from 
Ref. H3 that the d-symmetry HOMO would have lead to 
a high reactivity. However, the E^g HOMO is mainly of 
5s character and is not symmetry matched with the H 
a* orbital, hence it results in an extremely low reactiv- 
ity. Thus, the s-type HOMO of the Ag+PIC nanosys- 
tem effectively provides a shielding effect to protect the 
nanosystem from being attacked by the hydrogen atoms. 

Fig. @] shows the charge density difference of an Ag + 
PIC nanosystem 

Mr) = P A9+PIC (r) P Ag (v) J2 P«° m - 

a 

Here, p Ag+PIC (r) is the self-consistent electron density 
of the Ag + PIC nanosystem; p^ tom is the electron den- 
sity of the free atoms in the PIC- molecule; p Ag (r) is the 
density in reference system. The Ag sphere with radius 
15.12 a.u. is also schematically shown. The electron- 
density charge is close to z = 1.33 a.u. Fig. @] reveals 
that the electron density change Ap(r) is strongly lo- 
calized. We see also that the adsorbed PIC-molecule 
causes an accumulation of the reference charge on Ag- 
metal side due to the three carbon atoms with the coor- 
dinates (14.21, -13.71,1.33) a.u., (12.13,-12.03,0.99) a.u., 
(13.93,-16.37,1.31) a.u. and one hydrogen atom with the 
co-ordinates (15.25,-17.37,2.26) a.u. Based on the chem- 
ical bonding model of the carbon and hydrogen valences, 
the chemical structure of interaction between the Ag and 
(C,H) atoms of PIC-molecule has the form 

H-Ag-C = C-C-H 

The major part of electron density in C — C and C = N 
bonds is localized at the bottom of the band and have 
the sp and sp 2 type covalence of the chemical bonding. 
At the same time, the upper part of the occupied states 
contain mainly the bonding electrons which are concen- 
trated in Ag — H and Ag — C bonds and have (d, s)-(s,p) 
resonance type of bonding with the maximum displace- 
ment towards the Ag nanoparticle. 

V. CONCLUDING REMARKS 

Carrying out a numerical solution of the Dyson equa- 
tion (|A11I) in the Kohn-Sham density-functional method- 
ology we determined the electronic and chemical struc- 
ture of Ag+PIC nanosystems. The calculation of Agi25- 
particle embedded in a dielectric medium with the rela- 
tive permeability e = 61.1 (50% water solution of NaCl) 
shows that the amount of the electronic charge beyond 
the jellium edge (electronic "spill out") increases in the 
dielectric medium as compared to the vacuum and equals 
in this case 21.6 electrons. The resonant £2g-level [see 
FigE] arises at the energy 0.24eV under the LUMO-level. 
Both the Ag - resonances and e-i g states are found as so- 
lutions of Eq. (|A17l) with the full self-energy M(z). The 



electronic structure obtained for the Ag+PIC nanosys- 
tem has a terminal Ag — C bond. 

Three C atoms and one H atom take part in the ad- 
sorption of the PIC molecule on the Ag nanoparticle and 
accumulation Ap of charge takes place. The hybridized 
Ag — H and Ag — C bonds are of the (d, s)-(s,p) res- 
onance type with the maximal displacement of charge 
towards the Ag nanoparticle. 
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Appendix A: Green's-function technique to calculate 
the electronic properties of organic admolecule on 
metal cluster 

Our approach is based on the general concept formu- 
lated in Ref. H?]- It makes it possible to calculate the 
electronic structure of metallic substrates and to study 
the behavior of a compact cluster adsorbed on a nanopar- 
ticle. The principal theoretical tool is the scattering the- 
ory formalism, which considers perturbation of a metal 
substrate single electron potential by a spatially com- 
pact cluster, which is itself compact by virtue of screen- 
ing. It can therefore be treated as a localized scattering 
potential for Ag electrons of the substrate. In order to 
obtain the electronic structure of admolecule on an Ag 
nanoparticle a "matrix" scattering approach is adopted. 
Our objective is to construct the Green's function matrix 
Gfj,u(z) for the perturbed system, e.g., silver nanoparticle 
with organic PIC-molecule. 

The theoretical analysis providing us with a link be- 
tween the nanosystem and the related simple systems (Ag 
nanoparticle and PIC-molecule) is based on the Dyson 
equation ([7]). As is well known, the original problem 
of solving a linear differential equation can be mapped 
onto the solution of a matrix equation of infinite dimen- 
sion by expanding the wave functions in terms of a linear 
combination of properly chosen orthonormal functions, 
e.g. atomic orbitals (LCAO's) that are used here. In 
most cases of the LCAO cluster calculations Slater-type 
orbitals strongly facilitate the numerical calculation of 
overlap integrals. 

In practical calculations we construct the Green's func- 
tion matrix G fiv (z) © using a finite set of Nj, basis func- 
tions Xfj,( r )- The basis set only needs to cover the real 
space, within which Ap(r) is localized. We will denote 
this region as box A of volume Qa- In the present imple- 
mentation of our method the XnZm(r) are the Kohn-Sham 
orbitals 

Xnlm(r) =Vni(l r - R a \)Yl m (0 r -H a , <f>r-R a ) (Al) 

placed at appropriately chosen positions R a in the PIC 
assembled molecule. Here we choose the basis of atomic 
functions so that p = nml. ipnim(^) are the Kohn-Sham 
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radial wave functions obtained with the help of the po- 
tential AV(r), which in turn is calculated by means of the 
functions VnZro( r )- The selfconsistency iterations are re- 
peated until a desired convergence is achieved. Yi m (9, <p) 
are the spherical harmonics centered at R Q . 
Writing the Dyson equation 



G(z) = G°(z) + AG(z) 



(A2) 



we see that only the difference operator AG(z) (and the 
difference electron charge density Ap(r) related to it) 
need to be be actually calculated. This difference op- 
erator has the form 



AG(z) 



\-G°(z)-AV-(L-y) 



where 



AVau 



X*Jr)AV\p(r)] X v(r)dr 



and I is a unit matrix. Here the factor L 1 ensues from 
the assumption that the basis set x^( r ) was used in the 
Cholesky decomposition S = L • for the overlap matrix 



Suv 



in order to obtain the orthonormal basis. 

The density variation is calculated using the equation 

N b N b 

Ap(r) = Im 2 E A *V*m ^ ( r ) ( A3 ) 

/J, — 1 V — 1 



where 



CW = ((L- 1 ) t G°WL- 1 )^, 



A Pllv = ((L-^tApL" 1 ; 



and 



Ap 



AQ(z)dz. 



(A4) 



The lower integration limit e& is chosen in such a way 
as to include all the relevant Ag and molecule states, 
£homo is the HOMO energy. To compute the integral 
(|A4[) . we introduce a contour C in the complex plane z 
enclosing all the poles of the Green's function up to the 
highest occupied molecular orbitals energy in the charge 
density integration. 

Our computational scheme is based on the spherically 
symmetric Ag metallic nanoparticle Green's function 

G\vy-e) = Y j yUOA)G l {r,r'-,e)Y l * m {e',ct>'), (A5) 



decomposed in terms of the spherical harmonics^ Here 

Ri{r < ,e)Ni{r >1 e) 



r 2 Wi(e) 



(A6) 



with r < = min(r,r ),r> = max{r,r ). Ri(r < ,e) and 
Ni(r > ,e) are the regular and nonregular solutions of the 
radial Kohn-Sham equations with the potential V°(r) 
and energy e. Wi(s) is the Wronskian of the functions Ri 
and ./V;. We obtain Ri and iVj by integrating the radial 
Kohn-Sham equation, using the asymptotic behavior of 
these functions at the origin and infinity. As is generally 
known, the regular solution Ri behaves asymptotically 
as r l+1 at r — > 0. When r — > oo, the nonregular solu- 
tion Ni must be an outgoing wave for the continuum- 
energy region(e > 0), and it exponentially decreases for 
the bound-state region (e < 0). 

The Green's function of the nanoparticle is projected 
onto the localized basis 



(A7) 



I 111 



The following notations has been used above: 



(n\Q A \Mr<,e)) 



f X ;(r-R«)fll(r<,^r, 



(7V i (r > , e )|9 A |^) = <i/|e A |JVi(r>,e)>* 



(v\Q A \N l {r > ,s))= f xKr-R^Niir^dr, 



1, r £ Qa — volume of the 
box region A 

0, otherwise, 



Here we use the single-center method to compute all 
integrals in the Slater-type orbital basis^ The Slater 
orbital centered in a point, defined by its location vector 
R Q , is usually 



Xmm(r) = N n i | r — R Q |' exp(~C I r - R a |) 



Xyim(9r-R ,,4>r-R a 



(A8) 



The radial part of the Slater orbital is expanded over the 
Barnett- Coulson/Lowdin function (BCLF)^S 



r-Rn 



\n-l— 1 



exp(-C | r - R 



— ^(2\ + l)A n x -[ /2 (t,R a ,r)P x (^f) (A9) 



\fR a r ~^ 



R n r 
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where P n (z) are the Legendre polynomials of degree n 
and ^4.™~L 2 (C, J? a , r) are the BCLF's defined by the re- 
cursion 



A'l +1/2 ((,R a ,r) = -^A n x -[ /2 (£,R a ,r) 



with 



• / ^A+l/2(C) Rail 



l A+l/2 



(Cp<)k a+ i /2 (Cp>) 



where I\ +1 / 2 (z) and K A+1 / 2 (z) are the modified Bcssel 
functions of the first and second kind; the variables p< 
and p> stand for the min(R a , r) and max(R a , r) respec- 
tively. In the present implementation of our method the 
Kohn-Sham orbitals are placed at appropriately chosen 
positions R Q . 

Information on the Ag+PIC adsorption can be pro- 
duced from the Dyson equation ([7]) that describes the in- 
teraction between the free molecule and silver particle. In 
this Dyson equation the part of the unperturbed Green's 
function is played by the Augmented Green Function 



^Ag+PIC'( Z ) — 




(A10) 



In order to describe the molecule — silver nanoparti- 
cle interaction we use the approach known in the theory 
of transition metal impurities in semiconductor o 12 ^ 40 " — . 
The terms in the model Hamiltonian @ for Ag+PIC 
nanosystem are written as 



H 



PIC-Ag 



kin 



Mk 



h.c. 



and 



H 



PIC-Ag 



= V" W k k'c[c k ' 



kk 



Here, (di a ) and c\ (d\ ) denote the usual fermionic 
annihilation and creation operators, respectively, which 
are labeled by the indexes, k and i a , containing site, or- 
bital and spin degrees of freedom; Mki a is the s,d — i a - 
hybridization matrix element, and Wkk' is the matrix ele- 
ment for the Ag particle short-range potential scattering 
(the subscript "a" refers to the a-th adatom in the PIC- 
admolecule and i a stands for the corresponding electronic 
state and atom site). The Dyson equation reads 



Gi a iSz) = GZl(z)[l + M ia {z)G iaia {z)] (All) 



where 



GUM 



M ia (z) = M ia (z)Q- 1 {z) 



Off-diagonal elements even for the nearest neighbors are 
up to two orders of magnitude smaller than the diagonal 
ones, and are neglected. 

The PIC electron levels Si a are found self-consistently 
as solutions of the Kohn-Sham equations for the PIC re- 
lated orbitals in the nanosystem environment. The self- 
energy Mi a (z) contains two contributions. The term 



k 



(A12) 



describes the hybridization between the s, d-PIC-orbitals 
(if>i a ) and the Ag electrons (ipk) with the matrix element 

M kia = I Vfc(r)AF(r)^ (r-R a )dr. (A13) 

JQ.A 



The factor 



Q(z) = 1 - AV pot G° Ag (z) 



(A14) 



in Eq. (|A11[) describes the short-range potential scatter- 
ing, where 



AV pot = Y, 

kk' 



^(r)AV(r)W(r)dr. 



and 



G° Ag (z) = 



^(k\(z-H° L 



k 



LDA) 



l \k) = 



k 



£k 



(A15) 



(A16) 



is the single-site lattice Green's function for the elec- 
trons in the Ag host cluster described by the Hamiltonian 

n LDA- 

The Ag electron levels Ek are found self-consistently 
as solutions of the Kohn-Sham equations for Ag-related 
orbitals in the nanosystem environment. The Green's 
function (|All|) describes the hybridization between the 
PIC-electron orbitals and the electrons in the silver host 
cluster, where the Ag electrons are influenced by the po- 
tential scattering due to AV. If this scattering is strong 
enough, it results in splitting off of localized levels from 
the top of the band. This effect is also taken into account 
in Eq. (IA11[) : the positions of the corresponding levels 
before the hybridization are determined by the zeros of 
the function Q(z) in the energy gap A HOM °- LUMO of 
the Ag+PIC nanosystem, which arises in the electronic 
structure due to the potential scattering only. As a re- 
sult the equation for the deep level energy determined 
as a pole of the Green's function [Eq. (IA1 ![) ] within the 
framework of the LDA technique reads 



-e ia -AV±P A = M ia (z) 



(A17) 



It takes into account the resonance part of the scattering 
amplitude in the i a (PIC) channel and its mixing with the 
potential scattering states arising in the k (Ag) channel. 
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The imaginary part of the Green's function yields the 
spatial and energy electron distribution 



-ImY / ^-hi[(z-s ia -AV iaia )Q(e)-M ia (z)] (A18) 
7r as 



p(r,e) = lmG(r,r;e). 

Integrating over the energy we get the charge density 
distribution, whereas the local density of states, e.g., in 
the cell flA, reads 



ni oc (e) = / p(r,e)dr 
Jn A 

The change of the density of states is 

An(e) = Tr(G Ag+PIC ) ~ Tr(G^) 
or after straightforward calculations 
An(e) = 



The problem is treated self-consistently, starting with 
the trial set of LCAO Slater-type functions. The dif- 
ference potential in the zero approximation is a sum of 
the atomic potentials of the nanosystem. The selfcon- 
sistency procedure for AV^(r) is carried out in a mixed 
fashion. The first two iterations use the arithmetic av- 
erage scheme, which later on is effectively substituted 
by the Aitken scheme.— Just seven iterations produce a 
10~ 4 Ry selfconsistency. 
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